Ultra-high energy Neutrinos from Centaurus A and the Auger hot spot 
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The Pierre Auger collaboration has reported a correlation between Ultra-High Energy Cosmic 
Rays (UHECR) and nearby Active Galactic Nuclei (AGNs) within ~75 Mpc. Two of these events 
fall within 3 degrees from Centaurus A, the nearest AGN, clearly suggesting that this object is a 
strong UHECR emitter. Here we pursue this hypothesis and forecast the expected rate of ultra-high 
energy neutrinos in detectors like IceCube. In our baseline model we find a rate of ~ 0.4-0.6 yr~^ 
events above a threshold of 100 TeV, the uncertainty of which is mainly related to the poor knowledge 
of the physical parameters of the source and on the details of the model. This situation will improve 
with detailed high energy gamma ray measurements of Cen A by the upcoming GLAST satellite. 
This would make Cen A the first example where the potential of high energy mult i- messenger 
astronomy is finally realized. 

PACS numbers: 95.85.Ry, 96.50.S-, 98.54.Cm 
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The field of UHECR physics has probably taken a ma- 
jor step forward with the recent detection by the Pierre 
Auger Observatory of a spatial correlation between the 
highest energy cosmic ray events and nearby AGNs 
20 out of 27 events with energies above ~ 60 EeV corre- 
late with a nearby AGN within a radius of 3.1°. Further- 
more, 5 out of the 7 non-correlating events lie along the 
galactic plane where the AGN catalogues are incomplete 
and the largest magnetic deflections are expected. 

With such support for the hypothesis that AGNs are 
the main emitters of UHECRs it is timely to explore the 
consequences for other areas of high energy astrophysics. 
We will focus on the connection between UHECRs and 
neutrinos and investigate the possibility for detecting 
the neutrinos associated with the UHECR acceleration 
in AGNs. From the experimental point of view the de- 
tection of UHE diffuse and point source neutrinos in the 
km'^ IceCube detector is promising and the AMANDA 
collaboration has already reported interesting limits Q. 

Roughly 10 events are concentrated in the Centaurus 
direction, a region with a high density of AGNs, consti- 
tuting a hot spot in the Auger UHECR map. Our focus 
will be on Centaurus A (Cen A) as a case study (see also 
0). Two events fall near this galaxy, suggesting that it 
could be the first identified UHECR source. Indeed, Cen 
A, the nearest AGN at a distance of only ~ 4 Mpc 
has longbeen considered as a prime UHECR source can- 
didate [Tj . The problem of predicting the neutrino flux 
from Cen A is similar to the attempts to relate the ob- 
served UHECRs diffuse flux with a prediction, or at least 
with an upper bound, for the diffuse UHE neutrino flux 
[1, [l3] ■ We will employ basically the same approach 
for the case of Cen A with the help of the available data 
on its spectral distribution. Various models of neutrino 
emission from AGNs have been discussed in the past (see 
for example pd|). A recent update has been considered 
m We describe our model in more detail below. 



I. THE AUGER FLUX 

The expected number of events in the Auger array can 
be calculated starting from the total integrated exposure. 
The auger group reports S = 9000 km^ yr sr at present 
[l[. For point sources we need the exposure per stera- 
dian given by S/Jlgo where ilgo = tt sr is the Auger field 
of view corresponding to 60 degrees as maximum zenith 
angle. In addition the relative exposure lij{5) is required, 
weighting a source with declination 5 for the effective ob- 
servation time. oj{5) is parameterized according to [l3 | 
using 6 ~ —35° for the Auger declination and normaliz- 
ing to 1 the maximum. Assuming a power law shape for 
the energy spectrum F = Fo{E / Eq)^" we get 
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For the case of Cen A we have N = 2 events above a 
threshold Ec = 60 EeV with a source declination Sg — 
—47° and relative exposure uj{Ss) — 0.64 which gives 
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or E^F ~ 6 X 1022(£'/EeV)"-3 eVVm^ s. The uncer- 
tainty on the flux estimate is roughly a/2/2 ~ 70% 
from Poisson statistics. The intrinsic slope with just 
two events is very uncertain and we thus use a ~ 2.7 
as seen in the diffuse UHECR flux just before the GZK 
cutoff assuming it as generally representative of the typ- 
ical UHECR emitter. The uncertainty in the Cen A flux 
is therefore significant, but the situation is expected to 
improve as more statistics is collected by the Auger ar- 
ray allowing, in principle, to constrain the spectral in- 
dex directly from the data. Further, once the source is 
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clearly identified, also lower energy events can be used 
to reconstruct the spectrum, despite the larger magnetic 
deflections. We will see in the following, however, that 
the main source of uncertainty in the v flux is the AGN 
modeling rather than the UHECR flux uncertainty. An 
additional uncertainty is related to the possible system- 
atic error on the absolute energy scale of Auger of up to 
30%. An independent calibration is in principle possible 
exploiting the dip feature present in the UHECR spec- 
trum at ~ IQis eV dl. This gives = 80 EeV and a 
flux roughly a factor ~ 1.6 higher. Below we refer to this 
as the "dip" energy scale. 

Dropping the factor Q,qq/uj{5s) we can use Eq. ([2]) also 
to estimate the diffuse UHE flux. Using N , Ec, S from 
[111 we find F ~ 21(£;/EeV)"^-^l/km2 yr sr EeV or 
E^F ~ 2 X 1024(£'/EeV)°-3 eV^/ui^ s sr, in good agree- 
ment with the Auger estimate itself [l6j . 
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FIG. 1: Resulting neutron, u^+u^ and !/e+S'e spectra for an 
E~^''^ and 10^<-E/GeV<10^^ proton injection spectrum in- 
teracting on the photon field of Cen A simulated with SOPHIA 
in the optically thin source limit. 



II. AGN MODELING AND u FLUX 

To relate the expected neutrino flux to the observed 
CR flux several assumptions and an underlying model 
are unavoidably required. The usual scenario assumes 
that protons are shock-accelerated to ultra-high energies 
and then interact with ambient radiation or matter pro- 
ducing secondary neutrons and pious with an associated 
flux of gammas and neutrinos. If we assume that protons 
are magnetically conflned in the source so that only neu- 
trons can escape producing the observed flux of UHE- 
CRs, a direct link between neutrinos and UHECRs, or 
more generally an upper limit, is possible [^, [ll]. In 
the following we will use this hypothesis. The presence 
of particles accelerated up to 10^° eV indeed implies a 
magnetic field generally strong enough to confine the par- 
ticles themselves. Furthermore, even if the protons can 
finally diffuse out of the source, the acceleration region 
is generally expanding so that adiabatic losses limit the 
maximum attainable energies [25| . 

UHE protons can interact in the source both with mat- 
ter through p + p ^ p{n) + tts or with the local radiation 
field through p+7 p(ri)+7r''^'*'). Subsequent pion decay 
then produces photons and neutrinos. The coUisionless 
conditions required for efficient shock-acceleration imply 
a relatively low matter density so that at ultra-high en- 
ergies the dominant neutrino production channel is gen- 
erally the photo-hadronic one. We will thus focus on this 
process, using the Monte Carlo code SOPHIA [l^ to sim- 
ulate the interactions of protons in the Cen A radiation 
field and to normalize the relative yields and multiplic- 
ities of the secondary particles per interaction. These 
processes occur mainly close to threshold and produce 
only 1-2 pious per interaction, p-p interactions, instead, 
although disfavored, have quite higher ( > 10) pions mul- 
tiplicities [l3l and thus would also give higher neutrino 
multiplicities. Finally, we will neglect muons and pions 
synchrotron losses in the source MF that however are ex- 
pected to affect the i/ flux only at the highest energies 



(see [H [H for a throughout discussion on the role of 
MFs). 

The final neutrino yield has also a dependence on the 
Spectral Energy Distribution (SED) of the Cen A radi- 
ation field. In particular, if numerous enough ambient 
high-energy (x-ray/gamma) photons and thus high CM 
energies are available, in principle the multi-pion pro- 
duction channels can be activated and an higher neu- 
trino multiplicity per interaction can be achieved with 
respect to the expected low near-threshold yield. We 
consider especially i/ production in the nucleus for which 
detailed SED information is available. The jet case is 
more uncertain and will be discussed briefiy. A compila- 
tion of measurements of the Cen A nucleus is reported in 
[20t . The SED has an approximate double peak struc- 
ture with an infrared peak and a second soft-gamma 
peak. Upper limits in the TeV region at the level of 
few % of the Crab flux has been reported by the cur- 
rent generation of VHE cherenkov telescopes [U, . 
We adopt for the photon number density the crude ap- 
proximation of a broken power law n{e) cx e~i-9 over 
the range 0.001 eV < e < 100 MeV and n(e) cx const. 
for e < 0.001 eV, accurate enough to determine the ac- 
tual neutrino production regime. We show in Fig. [1] the 
resulting neutron and neutrino spectra for an E^^'"^ pro- 
ton injection spectrum interacting on the photon flcld of 
Cen A simulated with SOPHIA in the optically thin source 
limit, in which multiple proton-photon interactions have 
a negligible role. Indeed, it can be seen from the figure 
that the interactions occur near threshold for the whole 
relevant energy range, while only at very high, uninflu- 
ential, energies E,y > 10^° GeV we enter in the higher 
neutrino production regime. We therefore use the rel- 
evant average quantities valid in the E^ < 10^° GeV 
range, namely: ~ 0.1, (^„) ~ 0.5, as the fraction 

of proton energy transferred respectively to the neutrino 
{h'e + i^fi+ their antiparticles) and neutron component per 
interaction, and rj^n = {Eu) / {En) ~ 0.04 for the aver- 
age neutrino to neutron energy. The final neutron and 
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neutrino flux arc thus simply related via 
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where the factor {^^) / (^„) rj^n — 5 gives the average 
mean neutrino/neutron multiplicity. Notice that the 
multiplicity is > 3 implying that, even near threshold, 
more than one pion is on average produced per interac- 
tion. Our estimates are in fair agreement with ref. (isj to 
which we address the reader for a more detailed discus- 
sion of the various neutrino production regimes in photo- 
hadronic processes. 

If, on the other hand, the production site is located 
in the Cen A jet rather than in the nucleus, a softer, 
more x-ray populated, n(e) photon spectrum can be pos- 
sible [2^ and an higher neutrino multiplicity is achiev- 
able. However, in this case the site of acceleration would 
be the shock regions/hot spots in the jet, where a large 
scattering in the spectral indexes is observed [2^ so that 
a firm prediction is hard to establish. We also remark 
that Cen A is classified as a misaligned BL Lac with 
its jet pointing 20°-40°[25 away from our line of sight 
so that the relativistic boosting effect should not play a 
major role if the UHECRs come from the jet. 

To extrapolate the neutrino flux from the CR flux to 
PeV energies further modeling of the internal acceleration 
mechanism of the source is required. Extrapolating the 
spectrum to very low energies is clearly unrealistic 
and in fact several breaks in the slope of the energy spec- 
trum with subsequent steepening as the energy increases 
are predicted. We follow dOjUEl for the modehng of these 
breaks in the neutron spectrum and to relate it to the ob- 
served UHECRs spectrum. We thus consider a scenario 
in which an ambient proton spectrum cx Ep^''^ interacts 
with the low energy radiation field producing neutrons 
that, escaping from the source, decay into the observed 
CRs spectrum. The proton injection index 1.7 is in gen- 
eral agreement with the typical value « 2 expected from 
shock acceleration and it is chosen in such a way that the 
CR spectral index matches the value 2.7 at Ultra High 
energies (see below). Although the details arc generally 
quite model dependent two clear breaks are predicted in 
neutron/CR spectrum in the highest energy regime. In 
the first the spectrum steepens by one power when the 
pion production process becomes efficient while a sec- 
ond one power steepening is predicted when the source 
becomes optically thick to photo-hadronic interactions. 
The two breaks are generally close so we assume a single 
break at the energy £'{,. The resulting UHECR spectrum 
is then 
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There are thus three species in the model with differ- 
ent energy spectra: underlying, not directly observable 
protons, with an injection spectrum cx E~^''^ , photo- 
produced underlying neutrons, with a spectrum cx E~^''^ 
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FIG. 2: Final neutron, UHECRs and total neutrino spectrum 
in the model of this work, normalized to the Auger obser- 
vation at is ~ 60 EeV. Also shown is the fit to the MeV 
soft-gamma observations [2^ . [27j and the HESS upper limit 
around ~ 1 TeV [2l[. The Cen A UHE-7 curve, not shown 
for clarity, lies close to the neutrino curve. 



before the break and cx E^^-"^ after the break thus follow- 
ing the proton spectrum, and escaping "physical" neu- 
trons, with spectrum further showing one power steep- 
ening above the break, oc Eq^"^ . The escaping neutrons 
then decay back into protons far from the source consti- 
tuting the final CR spectrum whose UHE tail is observed 
in Auger. The various nuclear species spectra and the 
neutrino flux are shown in Fig. [21 The proton injection 
flux is not shown for clarity. Notice that the neutrino 
spectrum is in general supposed to follow the underlying, 
neutron spectrum via Eq. ([?]), i.e. unattenuated by pion 
losses and with a behavior cx above the break en- 

ergy. We will see however that the final expected rate of 
neutrinos is insensitive to the exact behavior of the spec- 
trum above the break. It also does not depend crucially 
on the slope below the break as long as it remains in the 
range between 1-2. The main parameter determining the 
neutrino rate is the actual break energy itself. To esti- 
mate this we use observations of the gamma spectrum 
from Cen A. 

The most important point to take into account to this 
aim is that gamma photons interact via pair production 
with the same low energy photon background relevant 
to the photo-hadronic interactions so that a break in 
the gamma spectrum can be related to a break in the 
UHECR spectrum. In particular, considering the ratio 
of the related cross sections and inelasticities it can be 
seen that d, Tn-yiE^) /T^-y{E^) ~ 4 x W'^En/E^, 
where Ti is the optical depth of the related process. 
When the source becomes optically thick to both pro- 
cesses (i.e. T„^ = T^^ ~ 1) we thus have the relation 
Eb-^ ~ 4 X 10~^ Ebn between the neutron and gamma 
break energies. Thus, differently from the determination 
of the neutrino multiplicity, the exact shape of the spec- 
trum is in this case crucial. Cen A observations in the 
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gamma band [2^, \ZA 123] show several breaks in the range 
100 keV-100 MeV with a photon spectral index a ~ 1.7 
for E < 200 keV and a ~ 3.0 for 10 < E < 100 MeV 
[2^. The exact energy of the highest break, however, 
is time dependent due to the intrinsic variability of the 
source. Further, the observations become photon limited 
above 200 MeV [l^ , making unclear if a further relevant 
spectral steepening is present above this energy. To be 
conservative we assume a gamma break at the highest de- 
tected energy bin found by EGRET at Et^ ~ 200 MeV, 
implying Ej^n ~ 10^ GeV. Anyway, given the importance 
of this parameter, that basically determines the normal- 
ization of the neutrino flux, in the following we also an- 
alyze the effect of a different choice for Ebn- We show 
in FigOour final neutrino, neutron and UHECR spectra 
together with a fit to the MeV soft-gamma observations 

m 



|27| and the HESS upper limit around 1 TeV [21 

We finally comment on the hadronic associated gamma 
flux expected to accompany the CR and neutrino fluxes 
(see also ref. [2^ for a specific analysis of the issue) . The 
and r]jn factors are indeed very similar to the neu- 
trino case so that the neutrino and gamma fluxes are pre- 
dicted to be very close in shape and normalization. How- 
ever, while neutrinos leave the sources just after produc- 
tion, gammas are subject to further processing through 
the development of an electro-magnetic cascade that de- 
pletes the high energy photon tail producing sub- TeV 
photons. Pair production in the low energy photon field 
of the source and electron synchrotron losses need thus to 
be taken into account for a prediction of the observable 
gamma flux. This, in turn, require further modeling of 
the source adding further uncertainties in the predictions. 
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III. EVENT RATE IN A NEUTRINO 
TELESCOPE 



Pion decay yields the flavor ratio v^'-v^'-Vt = 1:2:0. 
However, due to oscillations, we expect the ratio 
i/g : : I/,- = 1 : 1 : 1 at Earth. At energies Ei, > 100 TeV, 
neutrino telescopes are fully efficient both to tracks from 
charged current generated ^'s and to showers from i^e 
events, r leptons from interactions are expected to 
be detected both as showers near E = 100 TeV and as 
tracks in the higher energy range. For southern hemi- 
sphere detectors, like IceCubc, neutrino events from Cen 
A are down-going. In particular, Cen A, with a declina- 
tion 6 ~ —47° appears at a zenith angle of 43° in the 
IceCube field of view. 

Regarding the background, for E,y > 100 TcV the at- 
mospheric neutrino flux is negligible while the residual 
background of atmospheric muons has a quite steep spec- 
trum rapidly decreasing with energy. For simplicity we 
assume full efficiency for UHE neutrino detection, while a 
more careful evaluation would require a detailed Monte 
Carlo simulation. The final background is thus consti- 
tuted by the UHE diffuse i^-fiux itself whose relevance is 
linked to the detector angular resolution which at these 



TABLE I: Cen A event rate in IceCube for various break 
energies Eb and for the Auger and "dip" energy scales. 

energies is ^ few degrees 0]. We will limit ourselves 
to estimate the expected I'-flux while assessing the cor- 
responding statistical significance will eventually rely on 
the measured diffuse flux normalization. To calculate the 
event rate in a km^ detector like IceCube we assume for 
showers an effective volume of Vcs = 2 km'^. For track 
events the effective volume could be much higher due to 
the muon and tan range. However, for a zenith angle 
of 43° the available overburden is limited. Wc therefore 
conservatively use the same I4ff also for tracks. From the 
effective volume the expected event rate is 



r+00 

N = NapV,s / dE, F,{E,) 



(6) 



where Na is the avogadro number, p the density of the 
target material (ice in this case), Eth = 100 TeV is the 
threshold energy and a^g = 6.78 x 10-^^{E^/TcVf -^^^ 
cm^ is the CC cross section [sO] . 

Using the neutrino flux from the previous section, gives 
Af ~ 0.35 yr^^ or J\f ~ 0.56 yr^^ for the dip calibrated 
energy. Thus the conclusion is that IceCube should col- 
lect C'(few) events from Cen A in 5 years. This could 
be enough for a confldent detection of the source if the 
diffuse UHE i^-background is not too high. However, al- 
though the detection of Cen A may be challenging we can 
in principle extend the analysis to the whole Auger hot 
spot assuming that the related UHECRs emitters share 
the characteristics of Cen A. This results in a rate A/" ~ 2 
yr~^ in a region of radius about 10 degrees centered on 
Cen A which should be evident after a few years of obser- 
vations. As anticipated the estimated event rate does not 
dependent crucially on the neutrino flux slope while it is 
very sensitive to the break energy. Eh, which determines 
the relation between the neutrino normalization and the 
UHECR normalization as inferred from Auger. In Table 
U we show the scatter in the values of Af for different val- 
ues of Eb and for the Auger and "dip" energy scales. An 
order of magnitude variation in Af is in principle possible 
if the value of Ef, differs correspondingly by one order of 
magnitude. Clearly further observations in the gamma 
band would be desirable to have a more robust estimate 
of Eb- Fortunately the situation is expected to improve 
with the launch of the GLAST satellite that should pro- 
vide high quality data up to GeV energies and possibly 
beyond. Also, deeper observations from Cherenkov tele- 
scopes in the TeV range would contribute to improve the 
picture. 

Despite the uncertainties the prospect of neutrino 
detection from Cen A and its surroundings are quite 
promising with the exciting possibility to perform true 
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multi-messenger astronomy, observing for the first time 
a source in UHECRs, neutrinos and 7-rays. This would 
also allow for detailed studies of the source acceleration 
mechanism (3l| and, if flavor tagging can be achieved, 
neutrino exotic properties could be tested [l^, [13] . 

We conclude by commenting on a puzzling aspect of 
the Auger data [33|: Although many AGNs lie in the di- 
rection of the Virgo cluster, no events are detected. Al- 



though the statistics is low and this could be an exposure 
effect it is intriguing to notice that the issue can be set- 
tled by observations of the associated neutrino emission 
in IceCube. 
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